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Abstract 

In elasticity, the method of forces, wherein stress parameters are considered as the primary unknowns, is known as the 
Beltrami-Michell Formulation (BMF). The existing BMF can only solve stress boundary value problems; it cannot handle the 
more prevalent displacement or mixed boundary value problems of elasticity. Therefore, this formulation, which has restricted 
application, could not become a true alternative to the Navier’s displacement method, which can solve all three types of 
boundary value problems. The restrictions in the BMFhave been alleviated by augmenting the classical formulation with anovel 
set of conditions identified as the boundary compatibility conditions. This new method, which completes the classical force 
formulation, has been termed the Completed Beltrami-Michell Formulation (CBMF). The CBMF can solve general elasticity 
problems with stress, displacement, and mixed boundary conditions in terms of stresses as the p rimar y unknowns. The CBMF 
is derived from the stationary condition of the variational functional of the Integrated Force Method. In the CBMF, stresses for 
kinematically stable structures can be obtained without any reference to the displacements either in the field or on the boundary. 
This paper presents the CBMF and its derivation from the variational functional of the Integrated Force Method. Several 


examples are presented to demonstrate the applicability of the completed formulation for analyzing mixed boundary value 
problems under thennomechanical loads. Selected example problems include a cylindrical shell, wherein membrane and 
bending responses are coupled, and a composite circular plate. 
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Nomenclature 

strain energy in the Integrated Force Method (BFM) functional 
radius of shell 

complementary strain energy in the IFM functional 

modulus of elasticity 

material matrix 

plate or shell thickness 

plate or shell rigidity 

length of shell 

plate bending moments 

shell bending moment 

shell tangential force 

directional cosines of the outward normal 

components of surface tractions 

intensity of the distributed load 

radial coordinates 

boundary compatibility condition in terms of stresses 

traction conditions in terms of stresses 

temperature at the midsurface 

transverse displacement components 

prescribed boundary displacements 

potential of external loads in the IFM functional 
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Cartesian coordinates 

coordinates of shell contours 

coefficient of thermal expansion 

cylindrical shell parameter 

boundary of an elastic continuum 

temperature difference between inner and outer surfaces 

plane strain components 

Poisson's ratio 


p ’ p’ s 


G x?°y' T xy 

<Z> 




IFM variational functional for plane stress, plates, and shells, respectively 

plane stress components 

Airy's stress function for plane stress 

shell stress function 


particular integral for mechanical and thermal loads, respectively 
Q shell domain 


Introduction 

The method of forces, also known as the Beltrami-Michell Formulation (BMF), and its variant, the Airy’s stress function 
formulation, were the preferred tools of analysis in elasticity during the 1940’s and 1950’s. 1 * 2 In fact, solutions for many 
classical elasticity problems have been obtained via the method of forces. 1 “ 3 The method of forces, however, could not compete 
with the Navier’s displacement formulation, especially in analyzing plates and shells with displacement and mixed boundary 
conditions. Thus, the application of the method of forces diminished, and the displacement formulation gained popularity. The 
demise of the method of forces was not due to any intrinsic generic deficiency of the method but to the incompleteness of the 
formulation. Because a set of boundary equations was missing, the application of the classical BMF was restricted to solving 
only problems with stress boundary conditions. In other words, the Beltrami-Michell’ s force formulation can be used to solve 
stress boundary value problems, but it cannot solve the more prevalent displacement and mixed boundary value problems. The 
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missing set of equations, which completes the BMF, has been identified as the boundary compatibility conditions. At this time, 
these boundary compatibility conditions have been derived only from a variational formulation. Direct derivation of boundary 
compatibility conditions is not known, and this may be the primary reason why these equations were not formulated earlier. 
Augmentation of the classical BMF with these boundary compatibility conditions resulted in a novel force method — the 
Completed Beltrami-Michell Formulation (CBMF). The CBMF bestows equal emphasis on stress equilibrium and strain 
compatibility conditions. It is as universal as the Navier’s displacement formulation, solving all three classes of elasticity 
problems: stress, displacement, and mixed boundary value problems. Thus, the CBMF overcomes the deficiency of the classical 
BMF. The CBMF can provide solution to stresses without any reference to the displacements, either in the field or on the 
boundary, for kinematically stable structures. 

The primary purpose of the structural analysis is to determine the internal stress state in an elastic continuum. In the CBMF, 
stresses are obtained directly as a solution to a set of equations of this formulation. Displacements, if required, can be calculated 
from stresses using integration. In the Navier’s displacement method, displacements (whether required or not) must be 
generated first; then stresses are determined indirectly through differentiation. As a result, in the displacement method, stresses 
can become inaccurate, especially when approximate techniques are used. In the CBMF, problems with thermal and initial 
strains are handled directly by the compatibility formulation, whereas in the Navier’s displacement method, they have to be 
treated indirectly using the concept of equivalent loads. The development of the CBMF is further justified because all the 
solutions that have been obtained with the classical BMF have to be verified; that is, it must be determined whether the boundary 
compatibility conditions have been satisfied or not The noncompliance of boundary compatibility conditions for a classical 
elasticity solution is indicated in Ref. 4. 

The novel boundary compatibility conditions, the key ingredient in the CBMF, were accidentally derived and then 
identified during the formulation of the variational functional 5 of the Integrated Force Method (IFM) for the finite element 
discrete analysis. The IFM for the finite element analysis, which can be considered the discretized version of the CBMF, actually 
was formulated before the CBMF. Henceforth, in this paper, the force method for analyzing boundary value problems in 
elasticity and for analyzing plates and shells is called the Completed Beltrami-Michell Formulation (CBMF). The method of 
forces for the finite element numerical analysis is still referred to as the Integrated Force Method (IFM). 

The boundary compatibility conditions were reported earlier for two-dimensional elasticity problems, 5 and stress analyses 
using boundary compatibility conditions were published for rectangular 4 and circular plates 6 in flexure for mechanical loads. 
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This paper includes the formal presentation of the CBMF and its application to analyzing circular plates and circular cylindrical 
shells subjected to both mechanical and thermal loads . Specialized equations for plates and shells are derived from the stationary 
condition of the IFM variational functional, and several mixed boundary value problems are solved to demonstrate the capability 
of the formulation. The first problem is a circular plate made of two different materials and subjected to thermomechanical loads. 
The solution of the plate example demonstrates the application of the CBMF to problems with displacement and interface (or 
jump) boundary conditions. Cylindrical shells are analyzed next. The shell examples demonstrate the use of the CBMF when 
membrane and bending responses are coupled. In addition, this paper serves as an initial, yet an unified and systematic, attempt 
to bring back the method of forces for analyzing general elastic continua. It is anticipated that the development of the CBMF 
for various shell structures, wherein membrane and bending responses are coupled, may become a significant avenue for 
research. 


Completed Beltrami-Michell Formulation of Elasticity 

The basic concepts of the method of forces (the CBMF being its specialization for analyzing elastic continua) can be 
initiated from the stress-strain law, which is universal to all analysis formulations. The stress-strain law that links stresses { o} 
to strains [e] through a known material matrix [G] can be written as 

{<t} = [G]{ £ } (1) 

The stresses in Eq. (1) must satisfy the state of equilibrium, and the strains must satisfy compatibility conditions. In other 
words, the stresses in the method of forces can be determined from the stress strain law given by Eq. (1) and (I) the stress 
equilibrium equations and (II) the strain compatibility conditions. Displacements are not essential for the determination of 
stresses. 

A finite elastic continuum consists of a field and a boundary. Stresses and strains must satisfy equilibrium equations and 
compatibility conditions both in the field and on the boundary, respectively, as 

(la) Stress equilibrium equations in the field 

(lb) Stress equilibrium equations on the boundary (or traction conditions) 

(Ha) Strain compatibility conditions in the field 

(lib) Strain compatibility conditions on the boundary 
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In the method of forces, all equation sets (la, lb, Ila, and lib), including the compatibility conditions, are expressed in terms 
of stresses. The equation set of the classical BMF contained conditions (la, lb, and Da), but it missed the boundary compatibility 
conditions (lib). The CBMF utilizes all four conditions (la, lb. Da, and lib). 


Governing Equations for the Completed Beltrami-Michell Formulation 


Consider the CBMF equations in the following plane stress problem. For simplicity and clarity, homogeneous kinematic 
boundary conditions are considered, and initial deformations along with body forces are neglected. The derivation of the 
equations from the IFM variational functional for nonhomogeneous boundary conditions with body forces is given in Ref. 5 
and is not repeated here. However, a brief presentation of the IFM variational functional is provided in the appendix for quick 
reference. The equations, as obtained from the IFM functional, can be separated into five groups (la, lb, Ila, lib, and III) as 
follows: 


Group la: Equilibrium Equations in the Field 


da dr 


dx dy 


^ = 0 


dr d< 7 




dx dy 


- = 0 


Group lb: Boundary Equilibrium Equations (or Traction Conditions) 

^(cr) = <y x n x + r xy n y- p x =0 


S 2^) = r xy n x +ff y n y- P y =0 


(2a) 

(2b) 


(3a) 

(3b) 


where a x , a y , and are three components of the stress tensor; n x and n y are the direction cosines of the outward normal vector; 

and P x and P y are prescribed boundary tractions. In the field, the equilibrium equations are functionally indeterminate because 
three unknown stresses are expressed in terms of two (Group la) equations. 


Group Ila: Field Compatibility Condition 

The functional indeterminacy in the domain is alleviated through the field compatibility condition of St Venant, which 
can be written in terms of the strain components as 
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( 4 ) 



and in terms of the stresses as 


v2 K +<T y) = 0 


( 5 ) 


Equations (2), (3), and (5), in essence, represent the stress or the classical BMF in elasticity that was developed in 1900.* 


This formulation, which is incomplete, can only solve stress boundary value problems. 

Gr ou p Jib: Boundary Compatibility . Conditi on 

Three stresses on the boundary are expressed in terms of two traction equations, Eqs. (3a) and (3b), thus, there is one degree 
of functional indeterminacy. The field compatibility condition given in Eq. (5) alleviated functional indeterminacy in the field. 
However, because St Venant did not formulate the compatibility on the boundary, the stresses there remained indeterminate. 
The functional indeterminacy on the boundary, which made the Beltrami-Michell stress formulation incomplete, was alleviated 
by Patnaik 5 with the formulation of the boundary compatibility condition. This boundary condition, when expressed in terms 
of stresses for isotropic material, has the following form 


The set of three equations consisting of the traction conditions given in Eqs. (3) and the boundary compatibility condition 
given in Eq. (6) ensures stress functional determinacy on the boundary because three unknown stresses are expressed in terms 
of three equations. 

Equations (2), (3), (5), and (6) represent the CBMF, which ensures the functional determinacy of the stresses both in the 
field and on the boundary of an elastic continuum. The CBMF can solve a general elastic continuum with stress, displacement, 
or mixed boundary conditions. 

Group III: Continuity Con ditions (or Displacemen t Boundary Conditions) 

The stationary condition of the IFM functional, given by Eq. (48) in the appendix, also yields two displacement boundary 
conditions, and for the homogeneous case. 



( 6 ) 
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u = u = 0 {a) 


v = v = 0 (b) 


(7) 


where w and v are prescribed boundary displacements. In the CBMF, the displacement boundary conditions do not appear 
explicitly in the stress calculations, provided the structure is kinematically stable. The displacements, if required, can be 
calculated from stresses by integration using the kinematic boundary conditions. 4,6,8,9 

Completed Beltrami-Michell Formulation Solution Strategy for Composite Continuum 

The CBMF solution strategy for a composite elastic continuum with fields of Q\ and ^ 2 , and stress, displacement, and 
boundaries of r s , and T r , respectively (Fig. 1), are briefly described. 

Step 1 : Satisfy the field equilibrium and field compatibility conditions given by Eqs. (2) and (4), for both domains .Qj and 

(I 11 the displacement formulation, the Navier’s equations 3 have to be satisfied.) 

Step 2: Satisfy the traction boundary conditions given in Eqs. (3) and boundary compatibility condition given in Eq. (6) 
on contours r s and T uy respectively. (In the displacement formulation, equivalent traction conditions written in terms of 
displacements and displacement boundary conditions on contours r s and r ui respectively, have to be satisfied.) 

Step ID: On the interface boundary, F r , three conditions have to be satisfied: 
two residual equilibrium equations, 

2 1 1 (<7)-2*(a) = 0 (8a) 

2, n (cT)-2 2 a (<7) = 0 (8b) 

and one residual compatibility condition, 

c R i (<t)- c R ii (<t) = 0 (9) 

The functions 2(c) and <£( c ) were defined in Eqs. (3) and (6), and the superscripts I and II denote the domains C2\ and 
respectively. (In the displacement method at the interface boundary, two displacement and two traction continuity conditions 
have to be satisfied.) 

Step 4: Once the solution for stresses has been obtained, displacements, if required, can be calculated by integration. The 
evaluation of integration constants requires the kinematic boundary conditions. In the Navier’s formulation, the displacements 
must be calculated whether they are required or not Stresses are then calculated using the differentiation and the stress-strain 
law. 
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The composite structure can be solved by the CBMF or by the Navier’s displacement method. The problem, however, 
cannot be solved by the classical BMF because of the lack of boundary compatibility conditions for the boundary J F u and for 
the interface contour r r 


Properties of Compatibility Conditions 


Two properties of compatibility conditions for the case of a plane stress problem are given in this section. 

(1 ) The field compatibility condition, written in terms of displacement variables, u and v , becomes a trivial constraint, such 
as an identity [f(u y v)-f(u,v)] = 0, where / represents the field compatibility condition given by Eq. (4). The boundary 
compatibility condition given by Eq. (6), however, does not become a trivial equation when written in terms of displacements. 
In terms of displacements, the boundary compatibility condition given by Eq. (6) becomes 


d 2 V 11 

f d 2 u <?v 

n x + 

’ d 2 u 11 

fd 2 v d 2 u ) 

|<N 

1 

1 

'^dx 2 + dxdy j 

dxdy 2! 

{ dy 2 + dxdyj 


( 10 ) 


The nontrivial property of the boundary compatibility condition contradicts the popular belief that all compatibility 
conditions are automatically satisfied in the displacement method. 

(2) The field compatibility condition can be derived by eliminating the displacement components from the strain 
displacement relations. This logic as yet cannot be extended to derive the boundary compatibility condition. At present, the 
boundary compatibility conditions can be generated only from the IFM variational functional. This is , perhaps, a primary reason 
why the boundary compatibility conditions could not be formulated earlier. 


A pplications of the Completed Beltrami-Michell Formulation 

In this section, the CBMF is applied to the stress analysis of circular plates and circular cylindrical shells. Governing 
equations for both cases are derived from the stationary condition of the IFM functional. Several example problems are 
presented to demonstrate the CBMF solution procedure. 

Completed Beltrami-Michell Formulation for Bending of CircularJPlates 

The IFM variational functional for a circular plate subjected to mechanical and thermal loads is given in the appendix. Its 
stationary condition yields the following equations: 
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(a) Field equation of equilibrium: 


d 2 , x dM a> 

r +r « = 0 ai) 

dr v ' dr 

(b) Field compatibility condition: 

r|K - *f,) + d + - K) « < 12 > 

In Eqs. (11) and (12), M r and My are the radial and tangential moment, respectively; r is the radial coordinate; q is the 
intensity of the distributed load; h is the plate thickness; AT, a material constant, is defined as (£& 3 /12)(l — v 2 ), vis thePoisson’s 
ratio; a t is the thermal coefficient of the material; and At is the temperature difference between the upper and the lower surface 
of the plate. 

(c) Boundary conditions are specialized for various support conditions as follows: 

Simply supported contour 

M r = 0 (13) 

Clamped contour: 

(1 /K^M 9 - vM r ) + a t (At/h) = 0 (14) 

Note that the condition given in Eq. (13) represents the static boundary condition, whereas Eq. (14) represents the novel 
boundary compatibility condition. 

For the analysis of composite domains, transition (j um P) conditions on interfaces between regions made of different 
materials have to be established. These equilibrium and compatibility conditions at the interface follow: 



(15a) 


(15b) 


(15c) 


In Eqs. (15), superscripts I and II denote two regions of the composite plate made of different materials. The boundary 
compatibility condition given by Eq. (14) and the residual boundary compatibility condition at the interface given by Eq. (15c) 
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represent new equations for analyzing circular plates. These equations (which were missing from the classical BMF and are 
unique to the CBMF) make possible the solution of composite plates in terms of stress parameters only. 


In the CBMF, the thermal effects are accounted for on the right side of the compatibility conditions given by Eq. (12) and 
(14), whereas mechanical loads appear on the right side of the equilibrium equation given by Eq. (11). The Navier’s 
displacement method does not include the compatibility conditions in explicit terms, which is the rightful abode for thermal 
effects. The Navier’s formulation, however, accounts for the thermal effects in the equilibrium equations through the concept 
of work equivalent loads, which may introduce numerical errors when approximate solution techniques are used. 


The CBMF solution procedure is presented through the analysis of a composite plate (Fig. 2). The plate consists of two 
segments: an inner plate (JQ$ with radius a , material properties £; and V/, and thickness A,*; and an outer annular plate (£2 0 ) with 
inner radius a, outer radius b , material properties E 0 and v 0 , and thickness h Q . The inner plate (f2j) is subjected to a uniformly 
distributed mechanical load of intensity q, and the outer plate (Q 0 ) is exposed to uneven heating with the temperature difference 
At. The plate is clamped at the outer contour, given by r = b. This example illustrates the CBMF solution process for (a) using 
the boundary compatibility condition at a clamped contour, (b) analyzing composite domains by means of transition conditions, 
and (c) analyzing thermomechanical loads. 

Equations (11) and (12) are solved to obtain general expressions for the moments M r and for the regions Q and Q 0 , 
respectively: 


Example: Analysis of a Composite Circular Plate Subjected to Mechanical and Thermal Loads 



(16a) 



(16b) 



(17a) 



(17b) 
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where B\ % C\,D\ t £ 2 , C 2 , and Z> 2 are integration constants. 


The thermal load does not appear in the homogeneous solution given by Eqs. (17a) and (17b) because — (a At) = 0 . These 

dr 

six constants are calculated from the following six conditions: one boundary compatibility condition given in Eq. (14) at the 
outer contour (r = b ), three transition conditions given in Eqs. (15) at the interface (r = a ), and (3) two implicit conditions at 
the origin (r = 0). The implicit conditions require that the moments M r and My have the finite values at the origin. The solution 
is obtained for a specific composite plate with the inner plate (X2 Z ) made of aluminum and the outer plate (£2 0 ) made of steel. 
Numerical values for the material parameters in the domains and Q 0 , respectively, are taken as E z = 10.6 x 10 6 psi, v z = 0.33, 
= 12.6 x 10*Y°F, E 0 = 30.0 x 10^ psi, v 0 = 0.30, and a[°^ = 6.3 x 10~^/°F. The radii are a = 6 in. and b = 12 in., and the 
thicknesses are = 0.2 in. and h Q = 0.15 in.; the magnitude of the distributed load is q = 100 lb/in. 2 ; and the temperature 
difference is At = 50 °F. After the integration constants are determined, the final solution for M r and My is for the domain Q 
(0 <ir<a) 


and for the domain Q 0 (a<r<b) 


M'(r) = 844.05 - 20.81r 2 

(18a) 

A/^(r) = 844.05 - 12.44r 2 

(18b) 

M°{r) = 2046.63 - (5203.06 /r 2 )- 1 170 log r 

(19a) 

M°(r) = 2676.63 + (5203.06 /r 2 ) - 1 170 log r 

(19b) 


The displacements, if required, can be obtained by integrating the moment-curvative relations and by using displacement 
continuity conditions to evaluate constants of integration. Displacements for the domains and respectively, are given as 

w'ir) = 3.1209 - 0.0356^ + 0.1757xl0'V (20a) 

w°(r) = 5.3614 - 0.1344/- 2 - 0.7296 log r + 0.04417/- 2 log r (20b) 

The solution for M r My, and w (Eqs. (18) to (20)), respectively, has been verified from the corresponding solution with 
the Navier’s displacement method of analysis. 

Integrated Force Method Variational Formulation for Cylindrical Shells 

In this section, through an example of a circular cylindrical shell subjected to thermomechanical loads, the CBMF is 
extended to analyzing shell structures wherein membrane and bending responses are coupled. In CBMF, M x and force Ny are 
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the force unknowns for this problem. The variational functional for the problem is defined, and all equations of the CBMF, 
including the novel boundary compatibility conditions, for the cylindrical shell are obtained from its stationary condition. The 
EFM functional for the problem has the following form: 

n c s = A + B-W (21) 

where the strain energy A, the complementary energy B, and the work of external load W are given as 


-L 


M, 


d 2 w 


dx 


2 + ( 1 + v )«tT 




dQ 


(22a) 


_ f \JM X . . At\ dtyfHq, 


dQ 


(22b) 


W = f qw dQ 

JQ 


(22c) 


where w is the radial displacement; K = (Ef? /12)( 1 — v 2 ) is the rigidity; t a is the temperature at the midsurface of the shell; A t 
is the temperature difference between the inner and outer surfaces; and is the stress function. The stress function is defined 
through a procedure similar to that given by Washizu, 10 as 

M X = 'F (23a) 


d 2 * 


(23b) 


The variation flc has the following form: 


sn --L 


d M X 1 xr 

—r^ + ~ N <P +q 

dx* a v 


dQ8w 


-L 


M a . . At d 2 t 

K Eh dx 2 y ' ‘ h 1 dx 2 


dns* 


■\m x s(-—) + ^-Sw\ b +2 an 

( N <p 1 

-dr 

8 a 

i 

dN 9 

dtg) 

+ a t — 

SY 

[ x { dx) dx \ x 

a 

l\u ro J 

dx 

< Eh 

dx 

‘dx) 



( 24 ) 
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The stationary condition of the variational functional with respect to the displacement w and the stress function yields 


all the equations of the CBMF as follows: 
(a) Field equation of equilibrium: 


a v dx s 


(25) 


(b) Field compatibility condition: 


M a d 2 N & d\ 

K Eh dx 2 K 1 h ‘ dx 2 


(26) 


Contour terms in Eq. (24) yield boundary conditions. These are specialized for various support conditions as follows. 


Free Contour 

On a free contour, w * 0 and -7- w * 0, hence both the moment and its derivative must vanish: 

dx 


M x = 0 


sW- # 


(27a) 

(27b) 


Simply Supported Contour 

The rotation of the cross section is not prevented on a simply supported boundary. The condition S(dw/dr) * 0 results in 

M x = 0 (28a) 

Because w ± 0, the derivative -7"(aO = d*¥ / dx is not zero on a simply supported boundary, resulting in 

dx v 

(N 9 /Eh) + a t t 0 = 0 (28b) 

Equation (28b) represents a boundary compatibility condition for a simply supported boundary for the cylindrical shell. 


Clamped Contour 

For the clamped contour, both displacements and rotations are equal to zero, and on such a boundary two compatibility 
conditions must be satisfied. The boundary compatibility conditions, which are the coefficient of 5d}¥dx and d }¥ in the contour 
forms in Eq. (24), follow: 


(w t /ra)+<y, =0 


Eh dx 1 dx 


(29a) 

(29b) 
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Transition conditions, used to analyze composite shells, are derived similarly as those for circular plates. They consist of 
two residual equilibrium conditions 

M< I) -Af< II) =0 (30a) 

— + — M< n) = 0 (30b) 

dx dx 

and two residual boundary compatibility conditions given as 


-2~ +ah o 


n 

— +a f V 

E l f h ' ° 

1 1 J 


E 2 h 2 f ° 


= 0 


(30c) 


1 dUf j dt a 

— + a , - 2 - 

dx 1 dx 


dN i n) 




E 2 h 2 dx 


■ + a ; 


dx 


= 0 


(30d) 


The field equations given in Eqs. (25) and (26), together with appropriate boundary conditions represent the number of 
equations sufficient to solve the shell bending problem for stresses M x and The boundary compatibility conditions given 
in Eqs. (29) are derived for the first time. Two boundary compatibility conditions given in Eqs. (29) have to be imposed on a 
clamped boundary, and one compatibility condition given in Eq. (28b) must be satisfied on a simply supported boundary. Note 
that without the boundary compatibility conditions the solution ( M x and for the shell bending problems cannot be obtained 

for problems with either the displacement boundary conditions given in Eqs. (28) or the mixed boundary conditions given in 
Eqs. (29). The transition conditions given in Eqs. (30c) and (30d) enable the solution of composite shells by the CBMF. 

The field equations given by Eqs. (25) and (26) may be uncoupled to obtain the following alternative systems: 


d 4 M x 4 d 2 q Eh 

-f- + 4/TM = -- -y — T a 
dx dx a 


At d\ 


(1-t-v) — + a 


dx 1 


(31a) 




dx 2 


- + q 


(31b) 


or 


d N a , 

T- + 4p 4 N tn = 

dx 4 * 


Eh Eh 

q+ — a, 

Ka a 1 




h dx' 1 


dx H 


(32a) 
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(32b) 


„ ** d2N <P rr 

M. = t- + Ka, 


Eh dx 2 


/■> \ At 
(1+V) T 



where f$* = 3 (1- \^)/a 2 h 2 . Note that both Eqs. (31a) and (32a) are fourth-order equations, and either one can be selected for 
solution. Here, the moment equation (31a) is selected, and its general solution has the following form: 

M x = C { cosh fix + C 2 sinh fix + C 3 cos fix + C 4 sin j3x + (33) 

where Cj , C 2 , C 3 , and C 4 are the constants of integration, and and are particular integrals for distributed loads and 
temperature, respectively. The constants of integration are obtained by imposing appropriate boundary conditions. Once M x 
is known, can be calculated by back substitution from Eq. (3 lb). The solutions for two examples are provided to illustrate 
the CBMF solution process. The first example is a short cylindrical shell subjected to thermomechanical loads. The second 
example is a composite shell with clamped and simply supported boundary conditions. 


Example 1: Analysis of a Short Cylindrical Shell 

A simply supported cylindrical shell made of isotropic material with length L and radius a is shown in Fig. 3. The origin 
of the coordinate system is located at the centroid of the shell. The analysis is performed for two cases : (1 ) a uniformly distributed 
load and (2) uneven heating with At such that t 0 - 0. The material and the geometric parameters of the shell are such that the 
product pL<5; hence, it must be analyzed as a short shell. The general solution for each case is obtained by substituting the 
particular integrals into the general solution given by Eq. (33) and then by imposing the boundary conditions for simply 
supported contours at x = ±Z72 to evaluate the constants. 

Solution for the Mechanical Load 

For this case, = 0 and = 0. The solutions obtained for M x and after solving for the integration constants 
in Eq. (33) using simply supported boundary conditions at x = ±U2 have the following form: 

M x (x) = (q/2p 2 D){g 2 coshpx cos fix - sinh j 9x sin fix) (34a) 

%>(*) = i.9 a / D )[g 2 ($2 ” sinl1 P x sin P x ) + S\(&\ ~ cosh P x cos ( 34b ) 

where D = g 2 + g£ ,£i=cosh 2. cos X,g 2 = sinh A sin X; and X - PL/2. If required, the displacement w may be calculated using 
the stress-strain relations as 
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m>(x) = qa 2 /DEh j [g 2 (g 2 - sinh Px sin Px) + g l (g l - cosh (3x cos /}*)] 


(35) 


Solution for the Thermal Load 

For this case, particular integrals 'F^ = 0 and —K{\ + v)a t Mh, and the final expressions for the internal forces 

take the following form: 

M x (x) = (_«?<*>//>)[*,(*! - C osh px cos Px) + g 2 (g 2 - sinh Px cos /)*)] (36a) 

W = 2/3 2 /d) (g 2 cosh Px cos Px - g x sinh Px sin Px) (36b) 

If required, w can be calculated as 

w(x) = jEDh^{g 2 cosh Px cos 1 3x - g l sinh Px sin px) (37) 

The solution for this simple problem required a boundary compatibility condition (Eq. (28b)), even for a simply supported 
boundary, and hence the example could not have been solved by classical BMF. 

Example 2: Analysis of a Long Composite Shell 

A composite cylindrical shell of radius a and length 2 L is shown in Fig. 4. This shell is composed of two regions with 
different material and geometrical properties. Region I, bounded by contours 1-1 and 2-2, has material parameters and Vj, 
and thickness /ij ; and region II, bounded by contours 2-2 and 3-3, has material and geometric properties E 2 , v 2 , and h 2 . The shell 
is clamped along contour 1 -1 , and simply supported along the contour 3-3. Both regions are subjected to a uniformly distributed 
load of intensity q. Region I is also subjected to a temperature change of A t. 

Total solution for the composite shell is obtained by superposing the two component solutions. Each component solution 
involves four integration constants (C] to C 4 in Eq. (33)); hence, there are a total of eight unknowns for the composite shell. 
The eight constants of integration are evaluated from the following eight conditions: two boundary compatibility conditions 
for the clamped boundary 1-1 (see Fig. 4) given in Eqs. (29), two boundary conditions for the simple supported boundary 
3-3 given in Eqs. (28), and four transition conditions at the interface 2-2 given in Eqs. (30). 

For simplicity, the long shell condition (the products of P X L and P 2 L) is assumed for both components. Consequently, the 
response for the composite shell can be obtained by superposing effects from the three boundaries: that is, from the simply 
supported boundary, fixed boundary, and interface boundary as shown in Fig. 4. 
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Response From the Fixed Boundary ("Contour 1-1) 


The local coordinate system is defined such that the axis x\ is placed along the axis of the shell, with the origin in the plane 
defined by contour 1-1 (see Fig. 4). The solution process is similar to that presented for the short shell, given as 


M x (Xy) = (sin p x x x - cos fa) + «f^° 

(38a) 

N 9 (x x ) = aqle P ' x ' (cos jSjX, +sin P x x x )-1 


(38b) 


Using Eq. (38b), the expression for the displacement is obtained as 

w(*l ) = [~qa 2 /Eh^e ^ X ' ( cos P x x x + sin P x x { ) - lj (39) 

Response From the Interface (Contour 2-2) 

The expressions for M x> N^, and w, defined for regions I and II, respectively, are obtained as 


m ^ > (x 2 ) =e~^ 2 (A] cos p x x 2 + B x sin P x x 2 ) + 

(40a) 

N ( V(x 2 ) =2pfe~ fi ' x * (-5, cos p x x 2 + A, sin p x x 2 ) 

(40b) 

W (1) = {-\/aEh)2P x e (-fij cos P x x 2 +A X sin P x x 2 ) 

(40c) 

A/^^Xj) =e~^ lX ’ i (Aj cos P 2 x 2 + B 2 sin P 2 x 3 ) 

(41a) 

N ( v w (x 3 ) =2 p\e'^ (~B 2 cos P 2 x 3 +A 2 sin p 2 x 3 ) 

(41b) 

= (- l/aEh)2p 2 e (~B 2 cos P 2 x 3 + A 2 sin P 2 x 3 ) 

(41c) 


where Aj , B\, A 2 , and B 2 are the constants of integration, and x 2 and *3 are defined separately for each region (Fig. 4). Four 
constants of integration are calculated by imposing four transition conditions given in Eqs. (30) along the interface contour 
2-2. The four transition conditions yield the following four equations to compute four constants of integration: 


>• 

1 

II 

1 

Cl 

(42a) 

P 1 (b x -a x )+p 2 (b 2 -a 2 )=q 

(42b) 

2 ^ B x + 2 ^ 2 B 2 =— — q 

E\ E 2 \ 2 Efa Efy 

(42c) 

&(+*•>)* 

(42d) 


The solution of Eqs. (42) yields the four integration constants: 
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(43a) 


A - W%' 4,) /J!(-2^A -ift 2 -‘V|) + (9/2K>-*)(^? “ft 2 )] 

A “(‘/OJ^r’A’tft 2 + ‘ft 2 + 2i ftft) + (?/2Xl-‘)(‘ft 2 -ft 2 )] (43b) 

ft =(V%^' , tft?(fti 2 -tft 2 )4(,r/2ft)(l-t)(22ft 3 + ^^ + p 3)j (43c) 

ft = (1 - tftl)- (,/2ft Xl - b)(2ft 3 4 ft 2 ft 4 #$ )] (43d) 

where D =2fc^^ 2 (ft 2 + /)| ) + (ft 2 + AjSf) and A' = (fjAjJ/^Aj). The integration constants given in Eqs. (43) are intro- 
duccd into Eqs. (40) and (41) to obtain M x and for both regions I and II. Then, displacement vv, which can be calculated 
following the procedure given earlier, has the following form: 

Response From Simply Supported Conditions (Contour 3-3) 


For this case, a procedure similar to that presented for the clamped edge effects along contour 1-1 is followed. The 
coordinate axis X4 is defined as shown in Fig. 4. Contour 3-3 is simply supported, and the conditions given in Eqs. (28) are applied 
to obtain the expressions for the internal forces: 


and w is calculated as 


M x {x 4 ) = (q/2Piy^smP 2 x 4 

(44a) 

N <? (*4) = qa(e~ PlX * cos P 2 x 4 - 1) 

(44b) 

’(*4) = (** 2 M)(l “ e ~^ X * c°s 02*4 ) 

(45) 


As mentioned earlier, solution for any point is obtained by superposition of the expansions given by Eqs. (38) to (45). 
Analysis of the composite shell with simply supported, clamped, and interface boundaries can be obtained using CBMF. The 
problem, however, cannot be solved using classical BMF because of the following missing boundary compatibility conditions: 
Eqs. (29a) and (29b) for the clamped edge, Eq. (28b) for the simply supported boundary, and Eqs. (30c) and (30d) for the 
interface boundary. 


Conclusions 

The completed Beltrami-Michell Formulation (CBMF), wherein stresses are considered as the p rimar y variables, is 
obtained by augmenting the classical Beltrami-Michell Formulation (BMF) with novel boundary compatibility conditions. The 
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CBMF, which can analyze stress, displacement, and mixed boundary value problems in elasticity, alleviates the limitations of 
the classical formulation, which could analyze only stress boundary value problems. All equations of the CBMF for analyzing 
circular plates and cylindrical shells subjected to both mechanical and thermal loads have been derived from the stationary 
condition of the IFM variational functional. Transition conditions required for analyzing composite plates and shells made of 
different materials have been established. The CBMF has been used to solve several stress, displacement, and mixed boundary 
value problems in elasticity. In CBMF, displacements, if required, can be calculated from stresses by back substitutions. The 
CBMF is a true alternative to the Navier’s displacement formulation. 
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A ppendix — Variational Formulation for CBMF 


This appendix briefly describes the derivations of the CBMF governing equations for problems of two-dimensional 
elasticity and bending of circular plates. 

Governing Equations for Two-dimensional Elasticity 

The variational functional of the IFM for a two-dimensional domain 12 bounded by the contour /'has the following form 

n P /=A + B-W (46) 

where 



(47a) 


B 


= h\ 

Ja 


d 2 <t> v x -v°y ( d 2 <P <7 y -v^x d 2 <P 2(1 + v) 
dy 2 E dx 2 E dxdy E 


d£2 


(47b) 


W 


"*Ja (»«•+ VW-Jr -i ( F J* F ,*) d r 2 

1 1 2 


(47c) 


The variation of the functional with respect to u , v , and is given as 




do, fa. 




d, d 


5u + 


y J 


fa m do 


xy . _ y 


+ —^- + 5. 


dx dy 


dv 


dQ 


+ A f 

Eia 


p-h-v^)^K-va,)-2(l«)^ 


dns& 


~ */,. [K n x + X xy n y ~ P x) du+ (*xy n x + °y n y ~ 

+ A| r [(« + u)s[o x n z + T^n y ) + {v + + <r y n y )]dr 2 


iL fK- v sK|K- v<T -K- 2 ( 1+v ) 


dy' 


-fa n y + -dT n *) 


dr 2 84> 


(48) 
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In Eqs. (47) and (48), jTj and J 2 are portions of the boundary, where traction and kinematic boundary conditions, 
respectively, are imposed; B x and B y are components of body forces; an overbar denotes prescribed quantities; and <P is the 
Airy’s stress function. The stationary condition of the variational functional yields the governing equations of the CBMF given 
in Eqs. (2), (3), (5), and (6). 

Governing Equations for Bending of Circular Plates 

For a circular plate subjected to distributed loads q and the temperature change Ar, 


-*r 


M, 


d 2 w ' 

dr 2 


,, ( 1 dw\ 

+M *{-ri; J 


rdr 


(49a) 


B 


- h : 




M r ~vM v 

K 


At 

■+a , — 
r h 


d ( M m — vM r 
dr y \ K 


At 

+ a, — 

1 h 


rdr 


(49b) 


J. 


W = 2n 



a 


(49c) 


where r a and are the radial coordinates of the plate. 


The variation of the IFM functional with respect to variables w and *F can be written as 




n . dM m 

_( r « r )-_ + ^ 


drSw 


•/: 


K - VM <p d 2 ( M q> - vM r | a t dAl 


dr 


K h dr 


rdrS'F 




Af _ - vM At} 

JZ 1+a,— \S¥ 

K 1 h J 


(50) 


The stationary condition of the functional yields the CBMF equations for a circular plate given in Eqs. (1 1) to (14). 
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Fig. 1 Composite elastic continuum 
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Composite circular plate subjected to a uniform load, g, and temperature, At. 
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Fig. 3 Short, circular cylindrical shell. 
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Fig. 4 Long, composite, circular cylindrical shell. 
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